Imaging extrasolar planets by stellar halo suppression in 
separately-corrected color bands 
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ABSTRACT 

Extra-solar planets have not been imaged directly with existing ground or space telescopes 
because they are too faint to be seen against the halo of the nearby bright star. Most techniques 
being explored to suppress the halo are achromatic, with separate correction of diffraction and 
wavefront errors. Residual speckle structure may be subtracted by differencing images taken 
through narrowband filters, but photon noise remains and ultimately limits sensitivity. Here 
we describe two ways to take advantage of narrow bands to reduce speckle photon flux and to 
obtain better control of systematic errors. Multiple images are formed in separate color bands of 
5-10% bandwidth, and recorded by coronagraphic interferometers equipped with active control 
of wavefront phase and/or amplitude. In one method, a single deformable pupil mirror is used 
to actively correct both diffraction and wavefront components of the halo. This yields good 
diffraction suppression for complex pupil obscuration, with high throughput over half the focal 
plane. In a second method, the coronagraphic interferometer is used as a second stage after 
conventional apodization. The halo from uncontrollable residual errors in the pupil mask or 
wavefront is removed by destructive interference made directly at the detector focal plane with 
an "anti-halo" , synthesized by spatial light modulators in the reference arm of the interferometer. 
In this way very deep suppression may be achieved by control elements with greatly relaxed, and 
thus achievable, tolerances. In both examples, systematic errors are minimized because the planet 
imaging cameras themselves also provide the error sensing data. 



Subject headings: instrumentation: adaptive optics 
planetary systems 

1. INTRODUCTION 

Optical imaging of extrasolar planets, whether 
from space or the ground, requires improved meth- 
ods for halo suppression. These need to work over 
a broad wavelength band to minimize photon noise 
in both the wavefront sensor and planet images. 
Currently the best controlled starlight halos are 
obtained at ^1.65 (im wavelength from both space 
and ground, with speckle structure at ^2.5 x 10~ 4 
of the star intensity at ~0.5 arcsec radius (Schnei- 
der, 2002, Close et al, 2003). Subtraction methods 
such as rotation of the Hubble Space Telescope 
(HST) and spectral differencing for ground tele- 
scopes (Racine et al, 1999) yield a detection limit 
for faint companions at ^10~ 5 contrast. Order-of- 
magnitude improvements of these methods may be 



techniques: interferometric — stars: imaging — 

adequate to image massive, young giant planets in 
wide orbits, but older and less massive giants will 
generally be much fainter, even when heated by 
the star. Fainter still will be terrestrial planets, 
thus an Earth in a twin of the solar system at 10 
pc would appear in reflected light at ~2 x 10~ 10 of 
the sun at 0.1 arcsec maximum separation. Detec- 
tion at such faint contrast will likely require new 
and more powerful techniques, to directly reduce 
the photon flux (and noise) of the speckle halo, 
while maintaining exquisite control of systematic 
errors. 

Most proposed optical configurations for deep 
halo suppression are inherently achromatic, in- 
volving a combination of very efficient apodiza- 
tion or coronagraphy (Kasdin et al, 2003, Guyon 
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2003 and Gonsalves and Nisenson, 2003) with ac- 
tive correction of wavefront errors (Trauger et al., 
2003), and, in some concepts, also of pupil am- 
plitude (Littman et al., 2003). For imaging very 
close to the star, where terrestrial and old gi- 
ant planets are likely to be found, the halo com- 
ponent from optical errors is set by lower-order 
terms in the wavefront aberration, and these can 
be corrected by deformable mirrors of high accu- 
racy but modest spatial resolution (Malbet, Yu & 
Shao, 1995). However, if wavefront measurements 
are made across the pupil to determine even low- 
resolution correction coefficients, they would need 
to be both finely detailed and highly accurate, be- 
cause any under-sampled high-order structure will 
alias to low order. 

This difficulty can be overcome if wavefront and 
amplitude errors are derived from focal plane mea- 
surements. Thus, improved measurement accu- 
racy for phase and amplitude corrections to be 
applied at the pupil is obtained by phase diversity 
techniques applied to focal plane images (Malbet 
et al, 1995, Jefferies et al. 2002). This method was 
used to measure in detail the HST primary mirror 
aberrations (Roddier & Roddier, 1993). All focal 
plane techniques require the weak halo fluxes to be 
sensed in limited wavelength bands, and multiple 
channels will be needed to improve the sensitiv- 
ity set by photon noise. When rapid corrections 
are based on wavefront measurements made in the 
pupil plane, the limit to wavefront accuracy set 
by photon noise results in a halo with very weak 
speckles containing just a few photons each (An- 
gel, 1994). It can be shown that focal plane mea- 
surements of the weak halo speckles themselves 
can give the same accuracy and level of suppres- 
sion, provided the same total bandwidth is used 
(Angel, 2003). Focal plane measurements have 
the advantages that accurate calibration is not 
needed, because the measurement is null, and non- 
common-path errors are eliminated if the same 
detectors are used to provide the deep integrated 
planet images as well as instantaneous wavefront 
data. For our application in which the star im- 
age must be blocked, it is convenient to derive the 
errors intcrfcromctrically rather than by phase di- 
versity, obtaining the complex amplitude of the 
weak residual halo with the aid of a reference beam 
derived from the blocked starlight (Angel, 2003). 
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Fig. 1. — Concept for a narrow-band corona- 
graphic interferometer to sense and suppress the 
stellar halo, (a) A narrowband stellar image enters 
the system from above. The halo and planet image 
are reflected from a pierced mirror to be relayed 
via a pupil stop and beamsplitter onto the photon- 
counting detector arrays marked D and D^. The 
diffraction-limited core is passed through the mir- 
ror and a spatial filter to be used for phase ref- 
erence and/or destructive interference of the halo. 
Active or passive phase and amplitude control el- 
ements (e.g. a deformable mirror or spatial light 
modulator) may be placed in the halo light path 
at the reimaged pupil (Lyot stop), or in the refer- 
ence beam, (b) part of Airy pattern recorded in 
a lab demonstration, with the reference wavefront 
slightly tilted so that constant phase in the halo 
would appear as horizontal fringes. The fringe dis- 
placement of successive Airy rings reveals the half 
wave phase shift between them. 



2. HALO SUPPRESSION IN NARROW 
BANDS 

Since multiple wavelength channels are desir- 
able for accurate wavefront sensing, it makes sense 
to explore the potential for passive and active cor- 
rections within each limited bandwidth channel 
to directly suppress the halo. New modes then 
become possible, because separate correction of 
diffraction effects and wavefront errors is no longer 
required. To implement such schemes, an ini- 
tial separation of the telescope image into mul- 
tiple wavelength star images would be made with 
cascaded dichroic beamsplitters or a field grating 
(Bonnet & Courtes 1962). Labeyrie (2002) has 
pointed out one approach, in which a conventional 
coronagraph is followed by a second stage reimag- 
ing device. In it the halo speckle phases are ac- 
tively modified to be all the same, so at the center 
of an intermediate quasi-pupil their energy is con- 
centrated in a central spot that can be blocked. 
In our concept for suppression, a coronagraphic in- 
terferometer is used to both measure and suppress 
the halo in a servo loop, as shown schematically 
in figure 1. The active elements to control phase 
and/or amplitude are placed at the pupil image 
formed in the signal arm, or in the reference beam 
arm to modulate the reference image. One such 
instrument would be placed at each narrowband 
star image. In operation at the telescope, the cor- 
rections to be applied are determined from mea- 
surements of the complex amplitude at each focal 
plane image pixel, obtained from images recorded 
before and after a phase shifts of tt/2 in uniform 
amplitude reference beams. 

Multiple instruments, each with two-dimensional 
arrays of active elements, may at first seem a 
daunting proposition, but suitable small MEMS 
dcformable mirrors and liquid-crystal spatial 
light modulators for intensity and phase control 
(Liftman, 2003) are becoming available, as are 
CCDs with gain for use as fast, photon-counting 
arrays (Mackay, 2001). Wc consider below two 
new modes of operation. The first is a phase-only 
alternative to separate apodization and wavefront 
correction. In the second, conventional apodiza- 
tion and wavefront correction is followed by a 
second stage in which the residual halo is further 
reduced by destructive interference. 



3. METHOD I: PHASE CORRECTION 
OF BOTH DIFFRACTION HALO 
AND WAVEFRONT ERRORS 

In this method, the full halo (including diffrac- 
tion rings) is suppressed not by apodization or 
Lyot stop, but by modulation of phase across the 
pupil. The basic idea is to superpose Fourier phase 
components across the pupil that create the com- 
plex amplitude at each spot in the focal plane 
needed to null out the measured value. This is 
not possible over the whole focal plane, as can 
be seen from symmetry considerations. Thus the 
complex amplitude of the diffraction halo is Her- 
mitian, while from pupil phase modulation it is 
anti-Hermitian, so correction can only be made 
over half the focal plane. Also, the suppression 
is more efficient if the outer radius is limited, so 
we adopt a "D" -shaped nulled region in the focal 
plane. The required pupil phase shifts are found 
by iterative Fourier transform of the complex halo 
amplitudes measured within the D, and are of or- 
der 1 radian. We envisage that the phase pat- 
tern needed to suppress diffraction would be fig- 
ured permanently on the surface of a deformable 
mirror, and the active control used to clean up 
residual errors. Phase shifts made by displacement 
to correct diffraction can be accurate only at one 
wavelength A , and will be systematically off at 
neighbouring wavelengths. For reasonably small 
changes in wavelength, the relative phase error is 
given by &-<j)/4> = (A — A )/A , therefore there will 
be a chromatic increase in the halo intensity that 
varies as ((A — A )/A ) 2 times the original halo 
intensity. 

As an illustration, we show this method applied 
to the HST in figure 2. Obscuration of the pri- 
mary mirror by the secondary and supports results 
in a complex diffraction pattern. The optimized 
pupil phase modification yields about four orders 
of magnitude reduction for monochromatic light, 
with only modest loss of energy (the Strehl ratio 
is 43%). We have modeled the reduced halo sup- 
pression for finite bandwidth of 10%. The contrast 
at radius 3.5A/D (0.15 arcsec at A = 0.5/im) av- 
eraged over the band is 3 x 10~ 7 , still 3000 times 
below the uncorrected diffraction level. 

Such deep suppression and good throughput 
close to the star compares favorably with con- 
ventional coronagraphy. HST's NICMOS coron- 
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Fig. 2. — Correction of the HST diffraction halo by 
pupil phase adjustment, (a) The diffraction pat- 
tern at wavelength 0.5^m over 8 decades, given 
ideal figure correction, (b) Pupil phase change re- 
quired to correct for diffraction, optimized for field 
radius 2.5 - 35 X/D. The contour interval is 7r/4. 
(c) Monochromatic PSF after correction. Com- 
panions at 10" 6 star are revealed at radial separa- 
tions in 0.2 arcsec. (d) PSF at spectral bandwidth 
of 10%. 

agraph with a hard field stop of radius 2.1 X/D 
allows detection of companions as close as 2.6 - 
3.7 X/D, but with suppression of only 2.5 x 10 -4 
(Schneider 2002). A more aggressive coronagraph 
for HST, including optimized correction of low- 
order optical errors, could yield strong suppression 
of 3x 10 -8 , but at larger radius (7.5 X/D) and with 
only 18% planet transmission (Malbet et al. 1994, 
1995). A valuable aspect of the new method is 
that its deep suppression of diffraction and scat- 
tered halo can be maintained by adjustments with 
deformable mirrors alone, avoiding the reliance of 
conventional coronagraphs on precise manufacture 
and placement of a pupil mask. It will be partic- 
ularly powerful for complex or obstructed pupils 
for which apodization and Lyot stop methods are 
not very efficient (Sabatke et al, 2003). For exam- 
ple, it could be used to correct the halos caused by 
gaps and phase steps at the segment boundaries of 
the James Webb and Keck telescopes. 

4. METHOD II: ACTIVE HALO NULLING 

For deepest suppression over the full annular 
aperture, apodization would be combined with a 
second operational mode we now describe. The 
idea is to suppress the residual halo directly, by 



destructive interference in the focal plane. For 
this mode, the interferometer of figure 1 would be 
configured as a coronagraph, with halo suppres- 
sion achieved in the first instance by apodization 
preceding the interferometer. Based on a Fourier 
transform of the star halo's measured complex am- 
plitude, the deformable mirror at the Lyot stop 
within the coronagraphic interferometer would be 
adjusted to reduce the halo from phase errors. 
Some faint speckles will remain, because of lim- 
its to the manufacturing accuracy of the field and 
Lyot stops and in the setting of the phase cor- 
rection device. Their complex amplitudes would 
be re-measured, and the reference beam would be 
modulated to create a corresponding "anti-halo" 
with the same amplitudes but opposite phases. 
Destructive interference would then further reduce 
the halo. A planet image would be interpreted by 
the system as a speckle, but could not be sup- 
pressed, since it is incoherent with the starlight. 
Techniques for complex amplitude control with de- 
formable mirrors and spatial light modulators are 
described by Gonsalves, 1997, Roggemann & Lee, 
1997 and Liftman et al 2003. 

The additional reduction achievable in this 
way depends on bandwidth because of the ra- 
dial stretching of speckles with wavelength. Re- 
membering that phase changes by 7r for a ra- 
dius change of A/_D c ff (figure fb), wc find I/Iq 
(7r(?!)/2) 2 (AA/Ao) 2 , where AA is the full bandwidth 
and (j) is the field radius in units of A/D c ff • Asff is 
the effective aperture for resolution after apodiza- 
tion. 

The value of this second approach is that the ac- 
curacy now required for additional wavefront and 
amplitude control is much relaxed, because the 
complex amplitude to be adjusted is that of the 
weak interfering beam, not the full starlight beam. 
Thus, no matter how faint the starting halo, if we 
desire additional suppression by a factor of 100, 
the accuracy required in measuring and setting 
the anti-halo is only ~0.1 radians in phase and 
10% in amplitude. Such modest requirements are 
readily achievable with liquid crystal modulators 
of limited phase resolution. 

Figure 3 shows an example of suppression at the 
level required for detection of close-by Earth-like 
planets. The case modelled is a 4m space tele- 
scope, heavily apodized for < 10~ 12 suppression, 
but degraded by residual amplitude and phase er- 
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Fig. 3. — Halo suppression by destructive interfer- 
ence in the focal plane, (a) Intensity of a star halo 
(log 10 grayscale running from 1CP 6 -10~ 12 peak) 
calculated for an ideal, deep, radial Chebyshev 
apodization (IEEE 1979) of an unobstructed cir- 
cular aperture, adjusted for a halo at 10~ 12 level, 
with a companion of 10~ 8 star and 0.2 arcsec sep- 
aration. This mask has transmission of 7.5% and 
doubles the diffraction-limited image width. The 
scale is shown for D = 4 m (-D e ff ~ 2 m) and 
A = 500 nm. (b) The halo modelled for pre- 
dominantly low-order residual errors in the mask 
(power spectrum as k _5//3 and strength described 
in text), (c) The halo now reduced by destructive 
interference with an "anti-halo", modeled for 5% 
bandwidth as described in the text, (d) As (c), 
but including the effect of a spatial light modula- 
tor with discrete pixels subtending X/3D. 

rors that result in a halo of2xl0 _7 at0.2 arcsec. 
When the focal plane nulling is applied, cancel- 
lation to a level of 3 x 10~ 10 is obtained, even 
allowing for 5% bandwidth. To minimize pho- 
ton and speckle noise, measurements from 6 to 10 
bands would be combined. Optimum settings for 
the global and local correction elements would be 
worked out together. In this example, the required 
accuracy for the "anti-halo" setting is > 100 times 
more relaxed than that assumed for the initial di- 
rect pupil correction (0.2 nm rms in wavefront and 
0.25% in intensity at 0.25 m scale). 



5. DISCUSSION 

In general, a servo correction system will be 
necessary to repeatedly correct small distortions 
that develop in the wavefront. Only a few photons 
per speckle need be recorded over the full spec- 
tral band before a useful correction can be made. 
Repeated corrections made in this way will ensure 
that the halo flux is at the minimum possible level, 
and since each speckle is corrected before it has a 
chance to build up to more than a few photons, 
speckle noise in each cycle is held close to the level 
of halo photon noise. 

This approach is valid for both ground and 
space telescopes. The main differences are in the 
achievable suppression limit and in the servo cy- 
cle time (Angel, 2003). For ground telescopes, 
the wavefront is continuously perturbed by evolv- 
ing atmospheric turbulence. While some evolu- 
tion can be predicted and compensated, there will 
remain an unpredictable component that causes 
speckles at the few photon level to develop on a 
scale of a few tenths of a millisecond, depending on 
telescope aperture and seeing conditions. The cor- 
responding halo brightness for telescopes of 8 to 30 
m aperture at ~1 micron wavelength is likely to be 
^10~ 6 -10~ 7 of the star. For optical coronagraphs 
in space, the relative stability of the wavefront will 
permit much better halo suppression. But once 
the halo is at a level below that of the zodiacal 
background, there is little advantage in going fur- 
ther. For a telescope of a few meters aperture 
imaging at ~0.5 microns, this limit corresponds 
to ^10~ 9 of the star. In the photon noise limit, 
limiting sensitivity to exoplanets varies as (halo 
intensity) x -^/(cycle time), and is comparable for 
these ground and space telescope examples. 

The level of suppression needed for ground tele- 
scopes should be achievable by method I with a 
single very fast-acting deformablc mirror. The 
method has the advantage, mentioned above, that 
complex diffraction spikes from very large seg- 
mented mirrors can be suppressed with little loss 
of light. For the much higher suppression needed 
in space, the active halo nulling method will be 
preferred. 

This research is supported by AFOSR contract 
F49620-01-0383 and the NSF under grant AST- 
0138347. 
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